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We consider how trapped molecules can be sympathetically cooled by ultracold atoms. As a 
prototypical system, we study LiH molecules co-trapped with ultracold Li atoms. We calculate 
the elastic and inelastic collision cross sections of 7 LiH + 7 Li with the molecules initially in the 
ground state and in the first rotationally excited state. We then use these cross sections to simulate 
sympathetic cooling in a static electric trap, an ac electric trap, and a microwave trap. In the 
static trap we find that inelastic losses are too great for cooling to be feasible for this system. 
The ac and microwave traps confine ground-state molecules, and so inelastic losses are suppressed. 
However, collisions in the ac trap can take molecules from stable trajectories to unstable ones and 
so sympathetic cooling is accompanied by trap loss. In the microwave trap there are no such losses 
and sympathetic cooling should be possible. 

PACS numbers: 37.10.Mn, 34.50.-s, 34.50.Cx 



I. INTRODUCTION 

There has been rapid progress in the field of cold and 
ultracold molecular gases over the last decade, driven by 
a diverse range of applications in physics and chemistry 
Polar molecules are of particular interest because 
they interact strongly with applied electric fields, and 
interact with one another through dipole-dipole inter- 
actions that are long-range, anisotropic, and tuneable. 
These properties, along with the exceptional control that 
is possible at low temperatures over all the degrees of 
freedom, make an ultracold gas of polar molecules an 
ideal tool for simulating strongly interacting condensed- 
matter systems and the remarkable quantum phenomena 
they exhibit Q. In low-temperature molecular gases it 
becomes possible to control chemical reactions using elec- 
tric and magnetic fields and to study the role of quan- 
tum effects in determining chemical reactivity Cold 
molecules are also useful for testing fundamental symme- 
tries, for example by measuring the value of the electron's 
electric dipole moment [ij , searching for a time- variation 
of fundamental constants 043, or measuring parity vi- 
olation in nuclei [8] or in chiral molecules Q. For these 
applications, a great leap in sensitivity could be obtained 
by cooling the relevant molecules to low temperatures so 
that, for example, the experiment could be done in a trap 
or a fountain [lOj . 

The bialkali molecules can be produced at very low 
temperatures by binding together ultracold alkali atoms, 
by eitherphotoassociation [lj, [l2j or magnetoassocia- 
tion |13l - llq |. A few specific species of other molecules 
are amenable to direct laser cooling to ultralow tem- 
peratures (TtJ- A large variety of useful molecules can 
be produced with temperatures in the range 10 mK to 
1 K by decelerating supersonic beams [18l420l | or captur- 



ing the lowest-energy molecules formed in a cold buffer- 
gas source [2l|, [22J. For many applications it is desir- 
able to cool these molecules to lower temperatures, and 
this could be done by mixing the molecules with ultra- 
cold atoms and encouraging the two to thermalize. This 
sympathetic cooling method has not yet been demon- 
strated for neutral molecules, but is often used to cool 
neutral atoms (23l - [25| . atomic ions [26| and molecular 

ions mm. 

For sympathetic cooling to yield ultracold molecules, 
the rate of atom-molecule elastic collisions, which are re- 
sponsible for the cooling, must be sufficiently high that 
the molecules cool in the available time. In practice this 
requires that both atoms and molecules be trapped, so 
that they are held at high density and interact for a long 
time. The easiest way to trap molecules is in a static 
electric or magnetic trap. However, static traps can con- 
fine molecules only in weak-field-seeking states; since the 
ground state is always strong-field-seeking, inelastic col- 
lisions can eject molecules from the trap by de-exciting 
them to lower-lying strong-ficld-sccking states. These 
traps are therefore unsuitable for sympathetic cooling 
unless the ratio of the elastic to inelastic cross section 
happens to be particularly large. Inelastic losses can 
be avoided by trapping ground-state molecules, but such 
traps are more difficult to realize. 

In this paper, we consider the sympathetic cooling of 
LiH molecules with ultracold Li atoms. Due to its large 
dipole moment of 5.88 D, its low mass, and its simple 
structure, LiH is an attractive molecule for studying the 
physics of dipolar gases and the electric field control of 
collisions and chemical reactions. A supersonic beam of 
cold LiH molecules has been produced [2j| and deceler- 
ated to low speed using a Stark decelerator [3(J. Ultra- 
cold Li is likely to be a good coolant for LiH because 
the closely matching masses ensures that energy is trans- 
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ferred efficiently in an elastic collision. Also, the low mass 
of Li ensures that inelastic collisions with non-zero angu- 
lar momentum are suppressed by a centrifugal barrier, 
even at relatively high collision energies 31]. We have 
prepared a magneto-optical trap of 10 10 Li atoms for the 
purpose of sympathetic cooling. 

We begin by calculating the elastic and inelastic cross 
sections for LiH + Li collisions. Then we calculate the 
trajectories of a set of trapped LiH molecules that have 
occasional collisions with a co-trapped cloud of ultracold 
Li. Our aims are to calculate how the molecular temper- 
ature evolves with time, to investigate loss mechanisms 
in different kinds of traps, and to establish how the ul- 
tracold atoms should be distributed so that the cooling 
is most efficient. 



II. SCATTERING CALCULATIONS 

We have carried out quantum-mechanical scattering 
calculations on 7 Li+ 7 LiH collisions on the potential en- 
ergy surface of ref. HI. The calculations are carried out 
using the MOLSCAT program [H[. We use full close- 
coupling calculations for the energy range of importance 
for sympathetic cooling, up to collision energies of 1 K, 
and coupled states (CS) calculations over an extended 
range up to 100 cm -1 . The calculations are carried out 
treating LiH as a rigid rotor, with rotational constant 
bus/he = 7.5202 cm -1 . Because of the deep potential 
well (8743 cm -1 ) and strong anisotropy, a large rota- 
tional basis set is needed. The present calculations in- 
clude all functions with LiH rotational quantum number 
j up to j max = 37. The coupled equations are solved us- 
ing the hybrid log-derivative/ Airy propagator of Alexan- 
der and Manolopoulos [34| with the propagation contin- 
ued to 500 A. 

The collision calculations treat the Li atom as struc- 
tureless. This is justified because there are almost no 
terms in the collision Hamiltonian that can cause a 
change in the Li hyperfine state or magnetic projection 
quantum number [35l . |36| . The (very small) hyperfine 
structure of the LiH molecule is also neglected. 

The results of close-coupling calculations for LiH 
molecules initially in j = 1 are shown in Figure [TJ As 
expected from the Wigner threshold laws [37] , the elastic 
cross section becomes constant at very low energy (below 
about 1 mK) and the inelastic cross section is approxi- 
mately proportional to E" 1 in this region. Above about 
10 mK the ratio of elastic to inelastic cross sections sta- 
bilizes at a factor of 5 to 10. 

Close-coupling calculations are carried out for fixed 
values of the total angular momentum J, and the re- 
sulting partial wave contributions are summed to form 
cross sections. The results in Figure Q] include contribu- 
tions up to J max = 10, and Figure [2] shows the individual 
partial wave contributions for < J < 6. There is no 
significant resonance structure in the inelastic cross sec- 
tions for J < 6, although shape resonances appear in 
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FIG. 1: (color online) Elastic (red) and inelastic (blue) cross 
sections from close-coupling calculations on 7 Li- 7 LiH with 
7 LiH initially in j = 1. Also shown (black) are elastic cross 
sections for 7 LiH initially in j = 0. 



the elastic cross sections for J > 2 and are particularly 
prominent for J — 4 and 6. 

The cross sections obtained from quantum scattering 
calculations on an individual potential energy surface are 
in general quite sensitive to small potential scalings be- 
cause of variations in the scattering length a. However, 
this sensitivity is much smaller in Li+LiH because of the 
relatively low reduced mass. The low-energy limit of the 
elastic cross sections shown in Figure [T]may be compared 
with the value a = Ana 2 — 3200 A 2 obtained from the 
mean scattering length a = 16.01 A, as defined by Grib- 
akin and Flambaum |38j . 

The higher-energy results from CS calculations are 
shown in Fig.[3J It is interesting to compare the inelas- 
tic cross section with the Langevin limit, which assumes 
that all collisions that cross the centrifugal barrier lead to 
inelastic events. The Langevin limit is shown as a dashed 
line in Fig. [31 and it may be seen that the inelastic cross 
section remains below this limit even at collision energies 
around 100 cm -1 . 

We note that the reaction LiH + Li — s- Li2+H is highly 
endothermic and so cannot occur at the low collision en- 
ergies of interest here. 



III. SYMPATHETIC COOLING SIMULATIONS 

Using the cross sections calculated in Section II, we 
simulate the sympathetic cooling of LiH molecules co- 
trapped with ultracold Li atoms. We consider three types 
of trap. The first is a static electric trap for molecules in 
the weak-field-seeking state (j, m) = (1,0). Here, elastic 
collisions with the Li atoms cool the molecules, whereas 
an inelastic collision transfers a molecule to a lower-lying 
high-field-seeking state causing it to be lost from the trap. 
The other two traps we consider are a microwave trap and 
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FIG. 2: (color online) Partial-wave contributions to elastic 
(upper panel) and inelastic (lower panel) cross sections from 
close-coupling calculations. 
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FIG. 3: (color online) Elastic (red) and state-to-state inelastic 
cross sections from coupled states calculations on 7 Li- 7 LiH 
with 7 LiH initially in j = 1. The dashed black line shows the 
Langevin limit for the total inelastic cross section. 



then make a random choice between these two processes 
according to this probability. If the collision is inelas- 
tic, the molecule is lost from the trap. If the collision 
is elastic, the molecule's velocity is transformed by the 
collision into a new velocity using a hard-sphere colli- 
sion model. The velocity vector of the atom is selected 
at random from an isotropic Gaussian velocity distribu- 
tion whose width is fixed by the temperature of the atom 
cloud, chosen here to be 140 /iK. The atom and molecule 
velocities are transformed into the centre of momentum 
frame, where the molecular momenta before and after the 
collision, p and p', are related by 



an ac electric trap, both of which can trap ground-state 
molecules so that inelastic losses are avoided. We note 
that sympathetic cooling in an optical dipole trap has 
been studied previously [39(. 



p'=p-2(p-e)e, (1) 

where e is a unit vector along the line joining the centres 
of the spheres. It is given by 



A. Cooling in a static electric trap 

We first consider the sympathetic cooling of LiH 
molecules in the weak-field-seeking (j, m) = (1,0) state. 
The simulation starts with a large set of molecules with 
a velocity distribution that fills the trap. Later, when we 
consider ac and microwave traps, we will track individual 
molecular trajectories in these traps, but we do not need 
do this for the electrostatic trap. Our aim is only to cal- 
culate the fraction of all the molecules that cool to low 
temperature without being lost from the trap and this 
is determined entirely by the ratio of elastic to inelastic 
cross sections as a function of the collision energy. 

For each collision we find the kinetic energy in the cen- 
tre of mass frame, look up the relative probability for 
elastic and inelastic collisions as given in Figure [TJ and 



e = yi^|2P +b (2) 

where b is a vector that lies in the plane perpendicular to 
p and whose magnitude is the impact parameter of the 
collision normalized to the sum of the radii of the two 
spheres. For each collision, b is chosen at random from a 
uniform distribution subject to the constraints b • p = 
and |b| < 1. The new momentum of the molecule, p', 
is finally transformed back into the lab frame and this 
momentum is used in the next collision. 

Figure |4] gives the results of these simulations, show- 
ing how the fraction of molecules remaining in the trap 
and the temperature of their distribution depends on the 
number of collisions that have occurred. The molecules 
have an initial temperature near 100 mK, and for every 
10 collisions their temperature falls by about a factor of 
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FIG. 4: (Color online) Cooling and loss in a trap for weak- 
field-seeking molecules. The temperature of the molecules 
(blue circles) and the fraction remaining in the trap (red 
squares) are plotted against the number of collisions. 



10. After about 20 collisions, they have reached a tem- 
perature of lmK, but only 0.7% of them remain. After 
28 collisions, the molecules have thermalized to the tem- 
perature of the atoms, but now the fraction that remains 
is only 4 x 10 -4 . We see that the ratio of elastic to inelas- 
tic cross sections is too small in this case for sympathetic 
cooling in a static trap to be feasible. 

Our simulations use collision cross sections calculated 
in zero field, even though the molecules are electrostat- 
ically trapped. An electric field can have a large effect 
on atom-molecule collisions, as recently demonstrated for 
collisions between Rb and ND3 in an electrostatic trap 
[ioj . Here, it was found that the trapping field increases 
the inelastic cross section. If a similar effect occurred for 
Li-LiH, it would strengthen our conclusion that sympa- 
thetic cooling is not feasible in the static trap for this 
system. 



B. Cooling in an ac electric trap 

To eliminate trap loss due to inelastic collisions, it is 
desirable to trap the molecules in their ground state. The 
ground state of every molecule is strong- field-seeking, and 
strong-field-seeking molecules cannot be trapped using 
static fields. One solution is to use an ac electric trap 
where the molecules move on a saddle-shaped potential 
that focusses them towards the centre of the trap in one 
direction, but defocusses them in another direction. By 
alternating the focussing and defocussing directions at a 
suitable rate, molecules are confined near the trap centre. 
Such ac electric traps have already been used to trap 
polar molecules in strong-field-seeking states [4l| - |43| , and 
also to trap ground-state atoms |44J,|45|, and sympathetic 
cooling of molecules in ac traps has been proposed Q 

m. 

Motion in an ac trap consists of a small-amplitude mi- 
cromotion at the switching frequency of the trap, super- 
imposed on a larger-amplitude, lower-frequency macro- 
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FIG. 5: Phase-space acceptance of an ac trap in one dimen- 
sion for 4 phases of the switching cycle: (a) start of focussing 
period, (b) centre of focussing period, (c) start of defocussing 
period, (d) centre of defocussing period. The shaded areas in- 
dicate regions of phase space that are unstable in an idealised 
head-on collision. 



motion. The stability of the molecules derives from the 
micromotion, and a collision which interrupts the micro- 
motion may put the molecule onto an unstable trajec- 
tory. This means that a molecule initially confined in 
the trap may be ejected by a collision even though the 
collision reduces its energy. Figure [5] gives a simple pic- 
ture of how this can happen. In each dimension, the set 
of all the stable molecules forms an ellipse in phase space, 
and this ellipse evolves periodically with the phase of the 
switching cycle [42| . Figure 03a) shows the ellipse at the 
start of a focussing phase, showing that the positions and 
speeds of the molecules are positively correlated at this 
phase. To illustrate what may happen in a collision, con- 
sider the case where molecules collide head-on (b = 0) 
with stationary atoms of the same mass. In this idealised 
case, a collision reduces the speed of a molecule to zero, 
leaving its position unchanged, as indicated by the ar- 
row in Figure [SJa). The molecule will remain trapped 
only if it is still inside the ellipse, so the shaded regions 
of the ellipse are unstable against collisions. The same 
arguments apply at the start of the defocussing phase, as 
indicated in Figure 03c) . Half-way through the focussing 
and defocussing phases (Figure[SJb,d)), the molecules re- 
main inside the stable region when their speeds are re- 
duced to zero and so there will be no collisional loss at 
these phases. The phase-space plots show what happens 
in only one dimension. In a cylindrically symmetric ac 
trap, the same plots can be made for the radial and longi- 
tudinal directions separately, one being half a period out 
of phase with the other. Thus when the focussing phase 
begins in the radial direction, the defocussing phase is 
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beginning in the longitudinal direction, and a molecule 
must not be inside any of the shaded regions if it is to 
remain trapped after a collision. 

This picture is, of course, a highly simplified one. The 
collisions do not reduce the speed to zero, and they can 
couple energy from one direction to another, which tends 
to increase the opportunities for loss. Nevertheless, we 
expect the same general conclusions to hold - a large 
portion of the trap's stable phase-space volume becomes 
unstable when sympathetic cooling collisions are intro- 
duced, and collisions are more likely to result in loss at 
the start of a focus/defocus phase than half-way through. 

Turning now to a complete simulation, we consider LiH 
molecules in a cylindrical ac trap consisting of two ring 
electrodes and two cylindrically symmetric end caps, as 
used in refs. [Ill, 42]. The square of the electric field mag- 
nitude in this trap is well approximated by the expression 



E 2 (z,p) = EUl + 2a 3 
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where z is the characteristic size of the trap, Eq is the 
electric field magnitude at the trap centre, and 03 and a 5 
are the coefficients in a multipole expansion of the elec- 
trostatic potential. Considering the same trap as used in 
ref. [13], we set 03 = —1.29 and 0,5 = 0.63 for the longitu- 
dinal focussing phase, 0,3 = 1.29 and ag = 0.44 for the ra- 
dial focussing phase, z — 4.55 mm, and E = 50kV/cm. 
We switch the trap between the two configurations at a 
frequency of 5 kHz with a 50:50 duty cycle. 

An ensemble of initially warm molecules evolves within 
the trap, each molecule having occasional collisions with 
a distribution of 10 10 ultracold Li atoms. We calculate 
the trajectories of many molecules moving in the trap 
by solving the equations of motion numerically using a 
Runge-Kutta method with a fixed time step. The force 
acting on the molecules is F = — VW where W is the 
Stark shift. For the electric field magnitudes considered 
here, the Stark shift of ground-state LiH is small com- 
pared to the rotational spacing, and is given to a good 
approximation by second-order perturbation theory: 
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Here, /i e is the electric dipole moment of the molecule 
and &LiH is the rotational constant (in energy units). 

We suppose that the atoms are trapped independently 
from the molecules, for example in a magnetic trap, and 
we give them a spherically symmetric Gaussian spatial 
distribution with a 1/e half- width w a — 3 mm, and a 
Maxwcll-Boltzmann velocity distribution with a temper- 
ature T a = 50 fiK. 

We first simulate trajectories without any collisions, 
starting with an initial phase-space distribution that is 



larger than the trap acceptance, so as to obtain a set 
of molecules that, in the absence of collisions, survive 
in the trap for 10 s. This set of molecules defines the 
phase-space acceptance of the trap and is then used for 
the full simulation including the collisions. This ensures 
that molecules are lost from the trap only as a result of 
collisions. After each interval of time At, the calculation 
of the molecular trajectory is stopped and the probabil- 
ity, P, of the molecule having a collision during this time 
interval is calculated. The value of At is chosen such 
that P <C 1 and we take P = 11 a{K) v At, where n is the 
local density of atoms, v is the relative velocity of the 
LiH molecule and Li atom at the time of collision, and 
a(K) is the elastic collision cross section at collision en- 
ergy K . A random number, r, is chosen from a uniform 
distribution in the interval from to 1, and a collision 
occurs only if P > r. When a collision does occur, the 
velocity vector obtained from the molecular trajectory is 
transformed using the hard-sphere collision model out- 
lined above. The numerical integration of the trajectory 
then continues using this transformed velocity. Since the 
number of trapped atoms is many orders of magnitude 
larger than the number of trapped molecules, we assume 
that the atom distribution is unaffected by the presence 
of the molecules. We also neglect collisions between the 
molecules, since their density is so low. 

The dashed line in Figure ^i) shows the fraction of 
molecules that survive in the ac trap as a function of 
time. We see that most of the molecules are lost due to 
collisions and that this loss occurs on two separate time 
scales. During the first 1 s, 94% of all the molecules are 
lost from the trap. Between Is and 10s the number of 
trapped molecules continues to fall, so that after 10 s only 
1% remain in the trap. It is surprising to find that the loss 
continues at long times since we would expect there to 
be a small region close to the origin of phase space that 
is stable against collisions. It appears that even these 
molecules are eventually being destabilized by the colli- 
sions. To investigate why this happens, we repeated the 
simulation with the atom temperature reduced to zero. 
The result is shown by the dotted line in Fig. E^i). Here, 
the loss at early times is the same as before, but after a 
few seconds of cooling the fraction remaining in the trap 
stabilizes at around 5%. When the atoms have non-zero 
temperature, the collisions cause molecules near the ori- 
gin of phase-space to diffuse away from the origin, even- 
tually ending up on an unstable trajectory. The atoms 
cool the hotter molecules, but they also tend to heat the 
coolest ones, and even when the atom temperature is only 
50 fiK the heating results in significant additional losses 
from the trap on a 10 s timescale. 

Next, to shed some light on why there is so much col- 
lisional loss in the ac trap, we simplified the simulations 
by neglecting terms in the electric field beyond the sec- 
ond term in Eq. ([3J . In this harmonic approximation the 
phase-space acceptance of the trap is maximized, and al- 
though this field cannot be realized in practice [IH it 
is helpful to make this approximation since the dynam- 
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FIG. 6: (Color online) Simulated molecule loss in an ac 
trap, (i) Fraction of molecules surviving as a function of 
time. Dashed (purple) line: real trap, T a = 50 /xK. Dotted 
(blue) line: real trap, T a = 0. Solid (red) line: ideal trap, 
T a — 50 /iK. (ii) Number of collisions that result in molecule 
loss as a function of the phase of the switching cycle, in an 
ideal trap. The zero of phase corresponds to the centre of the 
radial focussing period and the labels (a-d) refer to the phases 
depicted by the same labels in Figure [5] 



ics in such an ideal ac trap are well understood. The 
higher-order terms complicate the dynamics by introduc- 
ing nonlinear forces into the trap and coupling the axial 
and radial motions, and this greatly reduces the trap ac- 
ceptance. Our simulations show that neglect of these 
higher-order terms increases the acceptance by a factor 
of 4 in both position and velocity. The solid line in Figure 
HJi) shows the fraction of molecules that survive in the 
ideal ac trap as a function of time. In this case almost all 
the loss occurs in the first 1 s of cooling. The loss occurs 
from the outer regions of the trap and there is a 'safe' 
region around the phase-space origin where a molecule 
will remain inside the trap's acceptance for all possible 
outcomes of a single collision. Once molecules have been 
cooled into this region there are no losses. The fraction 
of all the initial molecules that remain in the trap after 
3 s is 38%, and there are no further losses between 3 s and 
10 s. These results conform to the intuitive expectations 
obtained from our discussion of Figure [5] Comparing 
the results obtained for the ideal trap with those of the 
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FIG. 7: (Color online) Simulated time evolution of the ra- 
dial phase-space distributions of molecules in an ideal ac trap 
overlapped with a 50 /iK atom cloud with a width parameter 
of ra a =3 mm. Only the 1st and 2nd terms in Eq. ([3]) are 
included in the expression for the electric field. The phase of 
the switching cycle is the same in each plot and, to within 
one switching cycle, the times are (i) 0.002 s, (ii) 0.5 s, (iii) 3 s 
and (iv) 10 s. 



real trap, we see that it is primarily the large non-linear 
forces that make the ac trap an unsuitable environment 
for sympathetic cooling. 

Figure [6jii) shows how the number of collisions result- 
ing in trap loss depends on the phase of the switching 
cycle, in the ideal ac trap. As expected from the discus- 
sion of Figure El the simulations confirm that collisions 
occurring at the start of a focussing / defocussing period 
are far more likely to cause loss than those occurring 
half-way through these periods. Similar results are ob- 
tained when the higher-order terms are included, except 
that the modulation observed in the figure is then not so 
deep. 

Figure [7] shows how the radial phase-space distribution 
in the ideal ac trap evolves with time. At early times the 
molecules fill the available trap acceptance, but as time 
goes on they congregate near the origin of phase space. 
After 0.5 s [FigureJTJii)] the ellipse has become dense near 
the centre and sparse elsewhere. Most molecules have 
had one or more collisions by this time and these col- 
lisions tend to remove molecules from the outer regions 
of the distribution, either by cooling them towards the 
centre or, as discussed above, kicking them out of the 
trap. As time goes on, almost all the molecules in the 
outer regions of phase space disappear and the molecules 
that remain are cooled into a small region near the ori- 
gin. After 10 s [Figure Efiv)] this cold distribution has a 
full width at half maximum of 0.3 mm in radial position 
and 0.9 m/s in radial speed. The time evolution of the 
longitudinal phase-space distribution is similar. 
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C. Cooling in a microwave trap 

An alternative way to trap ground-state molecules is 
to use a microwave trap, as discussed in ref. 46:]. The 
ground-state molecules are attracted to the electric field 
maximum of the standing-wave microwave field inside 
a resonant cavity. The trap depth is particularly large 
when the detuning of the microwave frequency from the 
rotational transition frequency is small, although this 
places a stringent requirement that the microwave field 
be circularly polarised in order to avoid multi-photon 
excitation to rotationally excited states [46J. Collision- 
induced absorption of microwave photons may also oc- 
cur in the trap, and again this unwanted process is far 
more probable when the detuning is small [47J. Here, 
we consider a far-detuned microwave trap for ground- 
state LiH molecules, operating at a frequency of 15 GHz. 
Since this frequency is very small compared to the ro- 
tational frequency (2bun/h — 445 GHz), and since the 
Stark shift will also be small compared to 6lih for all at- 
tainable electric field strengths, the Stark shift is given 
to a good approximation by Eq. Q , where E 2 is now 
the time-averaged squared electric field. We take the mi- 
crowave field to be the fundamental Gaussian mode of 
a symmetrical Fabry-Perot cavity, having a beam waist 
of 15 mm and an rms electric field at the trap centre of 
E = 40kV/cm. This is the field produced by coupling 
2.6 kW of power into a cavity whose Q-factor is set by 
the reflectivity of room-temperature copper mirrors. The 
trap has a depth of 500 mK and the simulation begins 
with a trap whose phase-space acceptance is completely 
filled. We simulate individual molecular trajectories in 
the microwave trap with collisions modelled in exactly 
the same way as outlined above for the ac trap. We use 
the field-free elastic cross section shown in Fig.[T] since 
this is insensitive to the microwave field |47j . and we 
neglect inelastic relaxation between field-dressed states 
because, for our trap conditions, the rate of this inelas- 
tic process is expected to be very low compared to the 
elastic collision rate [47j |. 

Each time a molecule has a collision with an ultracold 
atom, its energy is reduced. Nevertheless, it is possible 
for a collision to transfer energy between axial and radial 
motions so that it has enough energy in one direction 
to leave the trap. By running simulations both with and 
without collisions, we find that there is no additional trap 
loss as a result of the collisions. This is because the axial 
and radial motions in the trap are weakly coupled, so 
that even in the absence of collisions a molecule whose 
energy is greater than the trap depth eventually leaves 
the trap. 

Figure [5] shows how the distribution of molecules in the 
trap evolves with time as they cool to the temperature 
of the atom cloud whose width is w a = 3 mm. Each plot 
shows the radial position and speed of each molecule in 
the trap. After 0.1s very little cooling has occurred and 
the molecules have the full range of speeds and positions 
that the trap can accept. After 0.5 s it is clear that the 
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FIG. 8: (Color online) Simulated time evolution of the radial 
phase-space distributions of molecules in the microwave trap 
overlapped with a 140 /iK atom cloud with a width of w a = 
3mm. The cooling times are (a) 0.1s, (b) 0.5s (c) 2s and (d) 
10 s. 



molecules are accumulating near the phase-space origin 
as expected. They have small speeds and are confined 
near the centre of the trap. As time goes on, the accumu- 
lation of cold molecules continues. After 2 s the majority 
of the molecules have been cooled into the small area of 
phase space near the origin, but some molecules remain 
distributed throughout the phase-space acceptance. The 
distribution has separated into two components, one cold 
and one hot. After 10 s, 90% of the molecules are in the 
cold component, the remaining hot molecules form a halo 
in phase space around the cold ones, and there is a re- 
gion in between where there are no molecules at all. The 
molecules that are slow to cool are the ones that initially 
have large angular momentum about the trap centre. In 
the absence of collisions, these molecules cannot reach 
the centre, and if they cannot reach the centre they are 
unlikely to have any collisions. The molecules in the halo 
in Figure H{d) have particularly large angular momenta 
and so spend all their time in the far wings of the atomic 
distribution; they are unlikely to have collisions even af- 
ter 10 s. 

The two-component speed distribution is even more 
evident when the atom cloud width is reduced to 1 mm. 
In this cold distribution develops rapidly at the 

centre of the trap but the rest of the trap phase-space 
acceptance is filled with hot molecules, apart from a thin 
empty region separating the hot and cold distributions. 
When the atom cloud width is w a = 5 mm the atom- 
molecule overlap is sufficient that the molecules form a 
single-component speed distribution. In this case it is 
possible to give a sensible measure of the temperature. 
The mean kinetic energy is not a good measure because 
a few remaining outliers with high kinetic energy have 
a disproportionate effect on the mean. Instead, we trim 
the distribution by removing the 5% that have the high- 
est kinetic energy, take the mean, and then divide by |/c 
to obtain a temperature. Figure H] shows the result. After 
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FIG. 9: Temperature versus time for molecules thermalizing 
with cold atoms in a microwave trap. The atom distribution 
has a temperature of T a = 140 fiK and a width of w a = 5 mm. 
To obtain the molecule temperature, we take the kinetic en- 
ergy distribution, remove high energy outliers, take the mean 
and divide by |/c. 



atom density is low. After 1 s, only 2% of the molecules 
are cold, but after 10 s, 96% of them are cold. As the 
atomic density is increased by reducing the size of the 
cloud, the cooling rate at early times increases. However, 
the number of cold molecules obtained after a long pe- 
riod of cooling is lower with these smaller atom clouds. 
For both the 1mm and 3 mm atom clouds, about 25% 
of the molecules are cold after 1 s. After 10 s, 90% of 
the molecules are cold when w a = 3 mm, but only 42% 
when w a = 1 mm. As discussed above, the molecules that 
fail to cool are those that have large angular momentum 
about the centre of the trap. These results show that the 
most suitable choice of atom cloud size depends on the 
trap lifetime. If the lifetime is long enough, it is best to 
use a large atom cloud to maximize the number of cold 
molecules obtained. If the lifetime is short, it is better 
to use a small atom cloud to maximize the cooling rate 
and then to remove the molecules that remain hot, for 
example by lowering the trap depth. 




Time (s) 

FIG. 10: (Color online) Fraction of cold molecules in the sim- 
ulated microwave trap as a function of time, for three different 
atom cloud sizes: 1mm (red solid line), 3mm (green dashed 
line) and 5mm (blue dotted line). A molecule is classified as 
cold if its kinetic energy is less than ffcT, where T = lmK. 



10 s the molecules have cooled from an initial tempera- 
ture of 100 mK to a final temperature of 200 fiK, close to 
the temperature of the atom cloud. The mean number of 
collisions per molecule required to reach this temperature 
is 30. As the molecules cool they move into the densest 
part of the atomic distribution and, as shown in Figure 
[H the collision cross section tends to increase. As a re- 
sult, the cooling rate tends to increase gradually between 
100 mK and 1 mK, despite the fact that the collision rate 
is proportional to the decreasing speed. 

Figure [10] shows the number of cold molecules as a 
function of time for three different atomic cloud sizes. 
The number of molecules is normalized to the total num- 
ber in the trap, and a molecule is taken to be cold once its 
kinetic energy is less than IfcT, where T = 1 mK. When 
w a = 5 mm, almost all the molecules in the trap cool 
to low temperature, but the cooling is slow because the 



IV. CONCLUSIONS 

Molecules are most easily trapped when they are pre- 
pared in weak-field-seeking states, but then sympathetic 
cooling is feasible only if the ratio of elastic to inelastic 
cross sections is high. In the Li+LiH system with LiH in 
its rotationally excited state {j = 1), we find this ratio to 
be approximately 5 at a collision energy of 100 mK, grad- 
ually falling as the collision energy decreases and reach- 
ing 1 at about 30 /iK. This ratio is too small for sympa- 
thetic cooling to be effective, since cooling from 100 mK 
to 100 [xK results in 4 orders of magnitude of trap loss. 
To avoid inelastic losses in this system, the molecules 
need to be trapped in the ground state. This can be 
done using an ac electric trap, but in this trap collisions 
can transfer stable molecules onto unstable trajectories. 
This occurs because stability of motion in an ac trap re- 
lies on a specific correlation between position and speed, 
and this necessary correlation tends to be upset by col- 
lisions. Our simulations suggest that the resulting trap 
losses are too great for sympathetic cooling to be feasible 
in a realistic ac trap. 

Alternatively, ground-state molecules can be trapped 
in the electric field maximum of a standing wave mi- 
crowave field formed inside a microwave cavity. The 
microwave trap appears to be suitable for sympathetic 
cooling. We find that both the cooling rate and the frac- 
tion of molecules that cool depend on the degree of over- 
lap between the atom and molecule distributions. When 
the atoms are compressed to a small volume, the cooling 
rate is high at the trap centre, but a large fraction of 
the molecules do not cool because their angular momen- 
tum prevents them from reaching the centre of the trap. 
When the atom cloud is larger, more of the molecules 
cool but the cooling is slower. The typical time required 
for a large fraction of the molecules to reach ultracold 
temperatures is a few seconds. For ground-state LiH, 
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in a good vacuum, the trap lifetime will be limited by 
black-body heating of the rotational motion [48]. This 
black-body-limited lifetime is 2.1s at room temperature, 
rising to 9.1 s at 77 K (49|. This suggests that liquid nitro- 
gen cooling of the microwave cavity may be necessary in 
order to cool a large fraction of the molecules in the time 
available. Cooling of the cavity mirrors would also allow 
for a factor of 10 increase in the cavity Q-factor, and a 
corresponding decrease in the power required to obtain 
the same trap depth. Note that the black-body heating 
rate is considerably slower for many other molecules of 
interest [48|, [4{| , and then trap lifetimes of several tens 
of seconds should be attainable under good vacuum con- 
ditions. 

Our use of a hard-sphere scattering model may under- 
estimate the degree of forward scattering, thereby overes- 
timating the cooling rate since scattering in the forward 
direction docs little to cool the molecules. We will inves- 
tigate this in future work. We have used an unchanging 
atomic distribution in our simulations, but it is clear that 
the atoms could be used more efficiently. It would be bet- 
ter to compress the atom cloud gradually so that the size 
of the atom distribution matches that of the molecules as 



they cool towards the trap centre. This optimizes the col- 
lision rate by optimizing the atom density and the over- 
lap between the two clouds at all times. Compressing 
the atom cloud will raise its temperature, but at high 
densities the atoms can be evaporatively cooled which 
will in turn sympathetically cool the molecules to even 
lower temperatures. We have focussed on LiH molecules 
sympathetically cooled with Li atoms, but we expect our 
general conclusions to apply to a wide range of other sys- 
tems. 
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